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The rigorous mathematical treatment for the propagation of a radio wave from a Hertz- 
dipole-source current-moment around a finitely conducting spherical earth surrounded by 
a concentric electron-ion plasma can be expressed as a series of zonal harmonics. Such a 
solution to the problem was obtained for the terrestrial sphere without a concentric plasma 
many years ago (1904-1915). However, the summation of the series, even at long wave- 
lengths or low frequencies, was considered to be impractical and the wi^ll-known and, indeed, 
rigorous Watson transformation was introduced (1918). 

The Watson transformation led to the development of elegant mathematical techniques 
both rigorous and approximate for the evaluation of the fields of radio waves in the vicinity 
of the earth. However, it does not necessarily follow that the Watson transformation 
is the only way to achieve numerical mastery of the problem. Indeed, it also does not 
follow that the Watson transformation is the most efficient approach to the rigorous form 
of the theory of propagation, especially at long wavelengths. 

This paper demonstrates that the field of the propagated long wavelength radio wave 
(frequencies less than approximately 50 kc/s) can indeed be evaluated by a summation of 
a series of zonal harmonics. Whereas the number of terms could become quite large (of 
the order of 10 k\a where a is the radius of the sphere and ki is the wave number of the me- 
dium between the concentric plasma and the earth), the speed with which these terms 
can be summed on a large-scale computer offsets the complications introduced by the Watson 
transformation as to the rigorous mathematical solution of the problem. 

The detailed structure of the field in the absence of a concentric plasma is characterized 
by the quite regular behavior of the ground wave as a function of distance. Indeed, the 
steady decrement of the ground-wave field is modified only near the antipode, where an 
interference pattern or standing w^ave as a function of distance is noted because of another 
wave's traveling around the sphere in the opposite direction. 

The introduction of the concentric electron-ion plasma shell traps the waves leaking 
into space, wliere reflection from the plasma builds up traveling waves in the direction of 
increased distance from the transmitter. Thus, the series of zonal harmomics comprises 
individual waves which are traveling in the radial direction with respect to the center of 
the sphere and standing in the direction of increased angular distance around the sphere. 
These waves, when summed, build up the wave progressing in the direction of increased 
angular distance. Under special circumstances, standing waves can be noted. This is 
especially obvious near the antipode of the transmitter. 

The results of the computations indicate that full rigor can be achieved with comparative 
ease at frequencies less than approximately 50 kc/s, leaving only the assumed model for 
the transmitter and the propagation medium and avoiding the complications of the Watson 
transformation. 

1. Introduction 

The treatment of the propagation of electromagnetic or sound waves diffracted by spherical 
objects as a series of zonal harmonics was introduced in such early theoretical investigations as 
those of Lord Rayleigh [1904] and Debye [1909]. The early research work on the radio prob- 
lem, such as that of Poincare [1904, 1910], Macdonald [1903, 1904, 1911, 1914], Nicholson 
[1910, 1911], and Love [1915], attempted to use the theory of zonal harmonics, either by a direct 
summation or by forming an integral from the series, to explain the ' 'remarkable observed facts 



1 This work was sponsored by the Rome Air Development Center, Griffiss Air Force Base, N.Y., mider Tasks 4 and 5of D.O. No. AF 
30(602)-2488. Techniques for summing zonal harmonics and evaluating Bessel functions for ELF/VLF wave propagation were developed for U. S, 
Navy, CNR, order No. NA-46-59. 
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concerning the propagation of radio waves around the earth's curvature.'' The series, as 
formuhited at this time, was quite difficult to sum because of slow convergence for the radio 
problem. Love, liowever, succeeded in summing a large number of terms, but his method was 
not considered to be completely satisfactory since it involved graphical construction and 
estimates. 

The results of Love, Poincare, and Nicholson did not explain the propagation around the 
earth's curvature, since it was correctly concluded for the assumed model that the diffracted 
field decayed exponentially, exp [-(3 {kia)^d], where /S is a constant, ki is the angular wave num- 
ber, ki^cc/c at a frequency y=aj/27r with c, the speed of light, a is the radius of the spherical 
earth, and 6 is the angular distance measured at the center of the spherical earth, figure 2, 
from the transmitter. 

Since the experimental evidence indicated that the field decayed approximately as 1-^d 
or at most 1/d, where d = ad is the distance along the surface of the earth, March [1912] refor- 
mulated the problem as an integral to obtain an answer more consistent with observed facts, 
and indeed found that the field decayed approximately as l/\^ sin d^ljd around the earth's, 
curvature. Although the work of March [1912] contained questionable approximations which 
affected the conclusions, a new approach to the problem had been introduced. March con- 
cluded that the waves were propagated without the exponential attenuation and further 
concluded that the ^'remarkable propagation of radio waves around the earth was explained." 
Apparently in answer to a criticism of March's [1912] paper by Poincare (Poincare [1904, 1910] 
previously concluded the wave was damped exponentially with distance), von Rybczynski 
[1913] extended March's paper and concluded that the exponential factor was neglected by 
March. He further concluded that this factor would become significant at great distances. 
However, von Rybczynski concluded that since the experimental results at that time were 
inconclusive (i.e., the results involved short distances) the significance of the exponential 
factor had not been established. 

In view of the discrepancies in the results of March and others, and apparently at the 
request of van der Pol, Watson [1918, 1919] reformulated the integrals used by March by 
means of the well-known Watson transformation and again reached the conclusion that the 
waves decayed exponentially around the earth's curvature, the experimental evidence not- 
withstanding. Watson [1919] further postulated a reflector of concentric plasma of high 
conductivity and reached the conclusion for this theoretical model that the waves did decay 
approximately as IJ-yjd around the earth's curvature as indicated by the observed facts, not for 
the reason set forth by March, but as a result of trapping of waves by the ionosphere. The 
comment of Watson in the 1919 paper ''a consequence of its presence (the ionosphere) is that it 
places grave obstacles in the way of communications with Mars or Venus, if the desirability of 
communicating with these planets should ever arise" is interesting from today's (1961) 
perspective. 

The Watson transformation seems to be the basis for most theoretical work in subsequent 
years (1919 to date), and Watson thus becomes the '^prophet" of a ''new rigor" in the treatment 
of the propagation of terrestrial radio waves. Thus, for example, the works of van der Pol 
and Bremmer [1937, 1938, 1949] use the Watson transformation for the radio problem. The 
treatment of radio wave propagation in the space between the ionosphere and the earth as a 
waveguide at very low frequencies [Wait, 1957, 1960] has also exploited the Watson transfor- 
mation. Although the technique of zonal harmonics has been used extensively in the literature 
on, for example, electromagnetic problems associated with cylinders [Wait, 1959], the appli- 
cation to the terrestrial sphere for the radio problem is confined to the work of Love [1915]. 
This is probably a consequence of the delusion that the slow convergence of the series of zonal 
harmonics for the radio-terrestrial sphere problem is an intractable situation. 

Whereas the Watson transformation has very great merit from a theoretical point of view 
and has been exploited extensively in certain theories of propagation, it is appropriate at this 
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time, especially at the low frequencies or long wavelengths,^ to re-examine the earlier treatment 
based on zonal harmonics wliich was apparently successful in yielding correct conclusions. 
Indeed, the procedure outlined in this paper, based on recurrence techniques for the calculation 
of Bessel functions [Goldstein and Thaler, 1959], can, as will be demonstrated herein, yield 
results quite reachly at frequencies less than approximately 50 kc/s. The formulation ])re- 
sented in this paper is completely rigorous and assumes only the models of a Hertzian dipole 
source current moment and propagation media consisting of a finitely conducting spherical 
earth with a specified concentric plasma shell. The results can also be extended to stratified 
spherical earth and indeed to a continuously stratified [Johler and Harper, 1962] concentric 
plasma, figure 4, as a model ionosphere. 

2. Theory of Propagation 

The electro dynamic fields E, volts/meter,^ and H ampere-turns/meter described by 
MaxwelPs equations, 

VX^+Mo-^^H-0 

VXH-eo^-E=J^ (1) 

where J is the conduction or convection current, amperes/square meter, can be treated as a 
continuous (^=— oo to +°°) time-harmonic wave, E=E (co, D), 11=11 (co, D), at a distance, D, 
from some reference and a frequency, /=co/27r, or, 

E--=\E\ exp [iwt-ikD], (2) 

where the positive time function is assumed, and the angular-wave number of the medium is 
represented by k = o)r}/c, where rj is the index of refraction of the medium. Maxwell's equations 
are therefore wirtten, 

IjlJcJi+VXE=0 

{eio: + a)E-VXH=0, (3) 

where e, a, and jjlo are the permittivity, farad/meter; the conductivity, mhos/meter; and tlie 
permeability, henry/meter, of the medium of propagation respectively. 

The splierical coordinate system, figure 1, with the unit vectors, r, d, cj), is introduced for 
convenience, since the earth and the ionosphere, for example, present approximately spherical 
boundaries to the propagation media. The origin of the coordinate system is located at the 
center of the sphere such as illustrated, figure 2. There are six components of the electrod}- 
namic fields in such a coordinate system, Er, Ee, E^,, Hr, He, and H^, which can be calculated 
by a differentiation process from a single [Hertz, 1889; Debye, 1909] vector, 11, provided the 
scalar IT satisfies the wave equation in the medium with wave number, /:, ( — A'^ = /Xoico(eico + o-)) 

(v'+k')n=o, (4) 

where n=IIa for constant a, and the operator V^ is defined by eq (55) or, 



2 Additional general references to the literature of long wavelength radio waves not otherwise cited in this paper but listed under 
"references" are Campbell [I960]; Jean et al. [1961]; Pierce [I960]; and Wait [I960]. 

3 Eationalized inks units are used in this paper. 
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ii0= exp io:t)\ t — ^ -^-^ 



£'0=[exp (io)t)] 



b dd 



t— 1 "I d r ()w~] 



He=[e.pii.t)]^^^^[rll- 



(5) 




-^=^0 



Figure 1. General coordinate si/stems. 




Figure 2. Spherical coordinate systems for model 
terrestrial sphere of radius j r=a. 




Figure 3. Spherical coordinates for model terrestrial 
sphere of radius r = a, surrounded by concentric 
plasma from r = c ^o r = co . 




Figure 4. Spherical coordinates, illustrating exten- 
sion to stratification of terrestrial model sphere and 
extension to a continuously stratified concentric 
plasma. 



740 



where 11'^ reefers to vertical magnetic source clipole or scatter field and 11^ refers to vertical 
electric source dipole or scatter field, provided the problem can be reduced to a two-dimensional 
form, figures 2, 8, and 4, or, 



dn 

d(f) 



= 0. 



(6) 



Tlie radio problem treated in this paper is one such two-dimensional problem. 

2.1. Primary Field 
The model for the source or transmitter is specified by the primary Hertz vector, 

where, for purposes of calculations, the source dipole current-moment, IqI ampere-meters, is 

47r _ 

taken as 1^1= — '^3.34(10 ^) ampere-meters. The radiation field of such a transmitter can be 

written approximately, 



£:, -10-^/0/-^ sin d exp [io:t-ikiD\ 
or, over an infinitely conducting plane of infinite extent, figure 1, x?/-plane: 

£',-£;,-2(10-^)/oZ ^exp [io^t-ikj)]. 



(8) 



(9) 



The radiation field is not necessarily the most important field at long wavelengths. Thus, 
unless the distance, D, is very great, the induction and electrostatic fields are quite important. 
It is therefore necessary to employ the exact expression for the assumed Hertz vector for the 
fields based on the Hertz vector (7), and eq (5), 



E, 



-1 r . ^d^n. „dn- 

~ -i—- — a Sin e ^r^ +COS 
rb sin 6 ]_ ^"' 



dff' 



be 



(10) 



where n= Ho, or for the particular case of transmitter and receiver located on the surface of 
the sphere, r=b=a, figures 2, 3, and 4, 



^ _ -1 fd'n cos e mi 



(11) 



Er=—^ exp [-\-io)t—ikiD] h — 



1 2 , cos d 



+ 



sm^ - 4 sm^ - 



kiG cos^ - . 

2 cos d 



Skia sin^ - 2 sin - 4:kia sin^ - 



The particular case, d=T, figures 2, 3, and 4, is of interest, 

£,=iexpM-2*<x][i-ig|^], 



(12) 



(13) 



since the transmitting and receiving dipoles are located end-to-end. The radiation field is 
therefore zero, but the induction and electrostatic fields remain finite and indeed not negligible 
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at extra low frequencies. However, at frequencies greater than, say, one kilocycle, 



^r^2a^ exp [-2i^ia](y>ikc/s). 



(14) 



Similarily, the horizontal magnetic, H^, field can be written, 



H^ 



^ 



fjLQicoa 



exp [io:t — ikiD] 



cos- 



2 sin- 



2kia sin - 



(r=b=a) 



(15) 



which for the particular case, 6=^Tr, 



H^=0(d=T), 



(16) 



2.2. Media of Propagation 

The media of propagation, figures 2, 3, and 4, are characterized by their electrical constants 
which can be expressed in concise form as the wave number, k. Thus, for air, with index of 
refraction, rj, or dielectric constant ei,=7}l, 



or 



/t — /ti — 2 ^1? 
C 



^1 -Vi, 

c 



(17) 



(18) 



(€i~ 1.000676 at the surface of the earth, ei = l in a vacuum, 771 '-^ 1.000338 at the surface of 
the earth, 771 = 1 in a vacuum), where c is the speed of light, C'^3(10^) m/sec. (c^2. 997925 (10^) 
m/sec in a vacuum.) 

The angular wave number of the ground, 



k=k,-='' 



^_ I . (7/Xo( 



(19) 



where a is the ground conductivity and €2 is the relative dielectric constant (relative to a vac- 
uum), where the permittivity, €=€?eo, eo=l/c^ mo? and Aio=47r (10~^) h/m, the permeability of 
space. 

Tiie angular wave number of the ionosphere is more intricate. 



k=k. 



Vo\e' 



(20) 



Thus, there are four distinct propagation components with the complex index of refraction 
77 og corresponding to ordinary wave (0), and extraordinary wave (e), upgoing wave (i), and 
downgoing wave (r). The index of refraction of the electron-ion plasma can be introduced 
with the plane coordinate system, figure 1. A simultaneous solution of Maxwell's eq (3) 
together with the Langevin equation of motion of the electron, 



m ^+mgV+^oe(VXHm) + eJ^=0, 



(21) 



where 7n is the mass of the electron, e is the electron charge, t is time, V is the vector velocity, 
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und H,n is the earth's static magnetic field vector, figure 1. Colhsions of electrons with ions 
can be treated either from the viewpoint of classical magneto-ionic theory where g^v and v 
is a constant average value collision frequency, or from the more accurate Bolt/juaini tlieor3^ 
of ionized gases in which electron collisions are proportional to energy, u, v = v{u) [Phelps, 
1960] where (j=g(o) + cx)ff), f/(co — co//), (/(co) [Johlor and Harper, 1962] for a gyrofreqiKMicy, o)h, 
or the parameter g is complex and frequency-dependent. 

Carrying out the simultaneous solution of the Langevin equations together with MaxwelPs 
equations results in a Booker [1939] t^^pe quartic in the parameter, f, assuming the convection 
currents, eq (3), 

J=-NeV, (22) 

or, 

ad'+as^+a2^'+a,^+ao=0, (23) 

where the index of refraction, rj, 

rj'^^'+^m'ct>,, (24) 

and (/)j, figure 1, is the angle of incich^nce on the lower ionosphere plasma, and the coefficients 
di, dz, (i2, Gb\, do iJ^i't^ given in the appencUx. The quartic (13) has four roots corresponding to four 
propagation components of tlie index of refraction (24) , tjq^I, In special case of a single concentric, 
uniform plasma illustrated, figure 3, only the upgoing indexes, 770, e, ai'e of concern. The down- 
going waves usually result from higher level reflections and correspond to the indexes T/o.e, figure 4. 
The validity of the Q-L approximation [Budden, 1951] at long wavelengths has been cliecked 
numerically [Johler and Walters, 1960], 



where 



:-exp[±^^,], (25) 

CO 



tan ft =^, (26) 

V 



2 

\[_ 



and — = ^ y \^U 



where coh and w.v are the gyro and plasma frequencies defined in the appendix and v is the 
constant classical colhsion frequency. The plus ( + ) sign is taken for the ordinary component, 
and the minus ( — ) sign is taken for the extraordinary component. 

2.3. Heavy Ions 

The effect of the heavy ions, i.e., the effect other than the collisions with electrons, cannot 
in general be neglected at extra low and ultra low radiofrequencies, since the gyrofrequency is 
approximately 50 c/s for these ions, according to Ratcliffe [1959]. ^.n approach to this problem 
has been introduced by Pfister [1955]. This approach can be developed as a modification of the 
calculation of the index of refraction, rj, of the plasma, but the formulation becomes quite 
intricate. 

The number of positive ions, A^+, the number of negative ions, A^_, in addition to tlie num- 
ber of electrons, N=Ne are to be considered. The convection current components are, figure 1 
eq (22) : 

J,= -NeV,= -e[NeV/-N+V^^+N_V.^]=J/ +J^ +J,~ 

Jy=-NeVy=-e[NeVy-N^V+y+N^V^y]=J;+Jy^+Jy- 

J,^-NeV,=^-e[Ny/-N^\\^~+N.VJ]=J,^+J,-^+J^~. (28) 
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The classical magneto-ionic collision frequency types, (v can be described in connection with 
certain constants, Ci, C2, C3, . . ., 

^ _ mo ^ rriQ ^ iuq ^ m_ 

'-'I — »'e,o j_ ' L'i = p+,o — — — — 5 C7 = i'-,o I — ' ^10 — ^0,- 



mo m+ + mo m- + mo moH-m_ 

^_ m_ m_ m+ ^ _ 7?2+ 

^2-»'<'.- ^ 1 ^ ' ^5-1'+ - —-!—-' C8-»'-.+ ——[-—' (^U-VO,^ 



m_ m+ + m_ m- + m+ mo-hm-i- 

n m+ ^ rUe ^ rrie ^ m<, 

me-\-m+ m+ + me m--\-me mo + me \^^t 

where Ve.^ is the electron-neutral molecule average collision frequency, Ve- is an electron-neg- 
ative ion average collision frequency, etc., where the average velocity of like particles is assumed 
to be zero, and m^, mo, m+, m_ are the masses of the neutral molecules, electrons, positive and 
negative ions respectively. The components of the Langevin eq (21) can then be written 
(assuming the Langevin equation can be applied to such ions) , 

1 1 lip ''"B 



II tp llvp 



flip 7/lg 



III A- III 4. 



III -L. Ill _(_ 

i<.V+,,+^E,+^[V+,,Hm.y-V+„H^J + {C,+C,+a)V+,-C,V,,-aV_.-aV,„=0 

III _L. Ill _1_ 

iwV_,y+^E„+^[V^,,H„,,-V_,MmJ + (.C,+ Cs+C,)V_,,-C,Vo.y-CsV+„-C,V,„=0 
ic.V.„+^E,+^[V_,Mn.y-V^„H^J + (C,+ Cs+C,)V_,-C,Vo,-C,V+.-aV,,,=0 

(^10+^11 + ^12) V^O, 2/ ClOy-,y Wl^+,2/ t/i2Ke,2/ = 

iC,,+Cn+Cn)Vo.-aoV.,-CnV+,-C,,V,,,=0, (30) 

where the components, figure 1, of E, and the velocity components of V are designated by the 
subscripts x, y, z of the plane local coordinate system at the ionosphere boundary. 
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The sinuiltaneous solution of these eqs (30) with MaxwelPs eqs (3) utihzhig eq (28) results 
in two coupled matrix equations: 



(In 

«31 

«'41 




an 







(l\2 







(X82 





«13 
^23 
^33 




^63 




^93 






au 









«^25 



0-45 









^36 
^56 






an 





a^i 


ail 

^97 





^28 








■ 



%9 




^79 





VJ%~ 




-K- 




J% 




a; 




J% 




E, 




Jt 




E. 


• 


Jt 


+ 


E, 




Jt 




E. 




J-. 




E, 




J-y 




Ey 




J^- 




E,. 



=0 (31) 



and 



hzi ^22 &23 
. ^31 ^32 ^33 , 



E.. 



E, 



+ 



^'\*Jy\^V \ ^ y) 
,C{Jl+Jt+J-^) 



-0, 



where (7= — l/^o;eo, and the a^' ^^^d 6^^- are defined in the appendix. 
The two eqs (31, 32) are equivalent to the single equation, 



^11 «12 ^13 «14 ai7 

0^21 a22 ^23 a25 ^28 

«31 «32 «^33 ^36 (^39 

a^l (^44 ^45 a4o ^47 ^ ^ 

a52 ^54 ^55 ^56 ^58 

^63 au G^65 «^66 «^69 

^71 ^74 ^77 a78 0^79 

^82 ass ^87 «8s ^89 

aga age 097 ags agg 



C00(700(700 
0(700C^OO(70 
00(700f700C^ 




1 



1 




1 








1 




1 




1 



hn bi2 bis 

621 U22 ^23 

O31 O32 O33 , 



■Jl 

d y 

J\ 

J\ 
7+ 

d y 

Jt 
Jl 

J y 

J^ 

E, 

Ey 

E, J 



=0. 



(32) 



(33) 



The 12 X 12 matrix, eq (33), and the vector submatrices can be divided into block matrices 
as indicated by the dashed line, 



■^9X9 
.^3X9 



^9X3 
V3X3_ 



L-^3X1. 



-0. 



(34) 
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Then M-J + N-K=0. 

and PJ+Q.K = 0. 

If M"' exists, multiply eq (35) by P-M~^ and subtract from eq (36). 
Q-K-P-M-^.N-K={Q-P-M-^-N)-K=Q). 



Then, 



(35) 
(36) 

(37) 



But, 



Q= 



K 


bu 


hn 


ft 22 


bn 


^32 



b. 



,imdK= 



Ey 

E. 



and P-M ^-N can be computed numerically [Johler and Harper, 1962; Johler and Walters, 
I960]. Thus (Q—P'M~^'N) can be written as a matrix involving the complex number, f, 
eqs (23, 24). {Q-P-M-''N)'K=0 has a nontrivial solution if, and only if, \Q-P-M-''N\ = Q. 
This yields a polynomial in f which can be solved to find f [Johler and Walters, I960]. These 
results depend on the existence of M~^ where Mis the original 9X9 matrix, eq (31). If E 
exists, the nonexistence of M~^ would imply that eq (31) has no unique solution. 
Writing, 



then eq (38) becomes. 



bn — ocn 


^2-«l2 


613 — «13 






621 — «21 b22 — Oi22 


b23 — Oi23 


= 


Q- 


bsi — asi ^32 — «32 


&33 — «33 






dix for aLGT) 




b^ — a^aT — oii2 


g=-(iL 

h = —a2s 


(l=—a^s 

i =aT^— I— a22 


.7 = — «31 

l^ar^+araL 




f+a 


b e^+d 




e 


f+/ 0^+h 




inj 


k 




l 



-PAI-'N 



(38) 



-«33, 



(39) 



t(l-ic) + t'{-cj- id-kg) + tVa+ U+ bgi-Jic-jd-kh) + f (6M 



-\-hjg-{-ekc~fcj—idj—kag)-\-{aJl + bhi-\-ekd—fjd—kha—lbe), (40) 



or, 



where, 



f+a 


b 


ct+d 


e 


r^+/ 


g^+h 


k+3 


k 


I 


03 = 


l — ar' 




^3 = 


dTasi—O-TCti 


i+kai^ 



=a4f'+a3f+a2f+aif+ao 



(41) 



a2=la+lj^baLaT—aT'—ai^a2i+ka22, 

ai = baTCX23 + ba^a^i — ekaT—fara^i —fa rccis+kaaL 

aQ=aj'l + ba22.a^x — ekai^—ja2,iaiz+kaa2z — lbe. 



(42) 



A preliminary numerical evaluation of the eqs (28-42) based on estimated approximate 
values of collision frequency and ion density indicated ion gyrofrequencies which could affect 
the index of refraction at long wavelengths (1 c/s to 100 c/s). 
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2.4. Boundary Conditions 

The conventional procedure^ for the sohition of electromagnetic-wave propagation problems 
requires the description of the pro[)agation in the various media assuming each of tliese to be of 
infinite extent in space. Thus, figures 2, 3, and 4, the electromagnetic wave is determined by 
the wave equation iu the Hertz vector and the angular-wave number, 

(v2+^?)n-() 

(V'+kl)U=0 

(43) 

where the arable numeral subscripts, 1, 2, 3, distinguish the three media of propagation described 
by the wave numbers ki, k2y k^. The boundaries r=a, c . . ., figures 2, 3, and 4, are then 
introduced by equating the tangential E and H fields and perhaps the normal H fields im- 
mediately on each side of each boundary such that tliese fields are continuous across tlie 
boundary. 

The boundary conditions for the particular example, figure 3, comprising a finitely-eoiuhict- 
ing sphere, a concentric plasma shell of infinite extent and uniform composition will now be 
considered in detail. The tangential components of E and H fields, eq (5) derived from appro- 
priate Hertz vectors, 111,9,3 for the various media are equated at the boundaries. The medium 
of wave number, k^ is a magneto-ionic medium involving ordinary and extraordinary wave 
numbers ks,o, ^3,e which can be treated as constants so that the wave eq (43) and correspou(hng 
Hertz vectors, 113,0, ^3,e, can be used (see appendix). The distinction between the two upgoing 
(i) propagation components and the downgoing (r) propagation components is not necessary 
in this case since the medium of wave nuuiber ^3 is of infinite extent. Downgoing waves arise 
from additional boundaries such as those illustrated, figure 4. The boundary conditions 
[Stratton, 1941] can be written by expressing the equality of the tangential components of the 
field on each side and in the immediate vicinity of the boundary, 

-^0,1^-^(^,3(0,6) 

E0,i = E0^^(^O^e) 
■tl e, I =^ -1^6,3(0,6) 

E4> .1 = E(j, ,2 
He 1^=^ He ^2 

Ee,i = Ee,2 (44) 

Straightforward substitution of eqs (5) into (44) would result in two sets of uncoupled equations 
for the electrical and magnetical solutions. Therefore, it is necessary to define the coupling 
in accordance with certain rules of the magneto-ionic theory. This can be accomplished with 
the Q—L approximation [Johler and Walters, 1960], 

Qoh^i^ n--^^^ (45) 

770 VO COS ^0 

^^-i, P.^^^, (46) 

Ve Ve COS de 

where sin </)i = r7 sin 6, atul Qo, Qe ^ii'c the ordinary and extraordinary components of —Qme and 
where the rectangular cartesian coorduiates are related locally at the model ionosphere boundary 
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to the spherical components, 



r\ ^ -^ 



Ez Er 



-sm 0^ 



or, 



Ey Ed 



• V 



(47) 



E^=QmeEd^ Ti-Ee 

Er=PmeEe^ 7 ' Eq 



The coupUng can be introduced into the boundary conditions by writing 

E^^i = QEe,z 



(48) 



(49) 



(50) 



at r=e. Equation (50), and the equated tangential fields, eq (44), can be written in terms of 

— > 
the Hertz vector 11, assuming a vertical electric dipole transmitter, r H.^ located in medium of 

wave number ^i. 



~ (rng+rnf)=^ {rut) 

nf=nf 
|;(rnr)=|;(rnr) 



(51) 



/:? (ng+no =Mo*"Q;^ ^ £ (rnr) 



' r dr ' 






(52) 



where in the Q-Jj approximation, 

ymg=±^ / J 2110 



f Pm 



— sm 4)1 



and p=sin^<^2 is the complex direction sine squared of propagation, and the minus sign (— ) is 
taken with the ordinary wave, k3,o, and the plus (+) sign is taken with the extraordinary wave 
k3,e. 

The quantity p, for the case of a plane wave incident upon the ionosphere, can be identi- 
fied with a local angle of incidence, 0^, at the model ionosphere boundary. 



p-^sin^c^^. 



However, in the case of the spherical waves incident on the ionosphere at r=c, 



-'*'=['+(fc)T' 



kV 



(53) 



(54) 
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a complex angle of incidence, which can be written explicitly for each spherical wave fn, ^n, as 
shown in the appendix. 

The solutioji of the wave eq (43) when dn_ 

\ V P ?1+-^. 4 r-- ^ ^l+^^n=o (55) 

r dr L drj rsmdddl_ ddj ^ ^ 

depends upon the separation of variables, 6 and r, 

U=J(r)F(d), (56) 

whereupon, 



^=\=- 






r 2^1=1^ r fi 



(57) 



Fid) sin ^ /(r) 

in which Legendrc differential equation is satisfied by the Legendre function, 

F{e)=P, (cos^), (58) 

where, 

Po (cos6>) = l, 

Pi (cos ^)==cos e, 

P2 (cos6')=J(3cos2 6'-l), 

P3 (cos e) = i(5 cos^ ^—3 cos 6>), 



and the eigenva ues, 

X=7i(n+1), 

n=0, 1,2, 3, ... , (59) 

or, 

r/2P (z) dP {?') 

(l-s^)^^^g^-23^^+7i(n+l)P„(2)=0, (60) 

where s=cos B. 

Also, eq (57) fStratton, 1941], since 



W^ldk^ ^J2=0, (61) 



where z=\^f(r), ^=kr, which is Bessel's equation. 

Thus, the solution of eq (55) can be written as a series of zonal harmonics consisting of 
terms of the form {Pni'n(kr)-\-yatn{kr))XPn{cos 6), where /?„ and 7„ are constants and 

fn{z)=j^JnH(^), (62) 



f.(s)=^f^i^ii(^), 



(63) 



Jn+112^^) and H^l\n2{z) are Bessel and Hankel functions (solutions of eq (61)) of order 71 + 1/2 
and argument and the Hankel function H^^^ is of the second kind fWatson, 1958]. The f-waves 
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Figure 5. Terrestrial sphere cross section with con- 
centric plasma illustrating various types of progress- 
ing spherical waves in the three media of propagation. 
Outgoing r-waves and ingoing lA-waves. 



CENTER OF SPHERE 



vanish at r= co , figure 5, and the ;/'-waves vanish at r=0. In particular, 
nS^iz: (2'^+l)f.(^i6)^.(M^.(cos e), (r<b) 

til/ o n=0 



=pzi^{2n+l)Ukir)Ulcib)Pn (cos 6), (r>6) 

til/ U 71 = 

nt=^fj {2n+l)[b%Uhr)+ct,Mhr)]Pu (cos 6} 
^i^iij^i'^n+inb^rnM+c'^UMWn (cos e) 

tilt n=0 

ni=T;^ t: {2n+l)at,Uhr)Pn{cos 6), 
^2=J;t t: {2n+l)a^Uhr)P„{cos 6), 
ni=7iT 2 {2n+l)diUhr)Pn(cos 6), 

It^l O ji=0 

^^=]^ § ('^n+l)d^Uhr)Pn{cos e). 



(64) 
(65) 
(66) 
(67) 
(68) 
(69) 
(70) 



The subscripts, 1, 2, 3, refer to the Hertz vector fields in the various media, figure 3, and 
a%, 6^, etc., are constants determined by the boundary conditions, eq (71). 

The boundary conditions, eqs (51, 52), can be written in concise form as the matrix equation: 



"ai.i ai,2 0^1,3 

^2.1 C^2,2 ^2.3 



^3,2 «3,3 

^4,2 %3 
















«5.4 








^6, 5 ^6, I 










■ 



«3,8 
«^4,8 



%7 

«6. 7 





^7,5 (17, fi ^7,7 

^8,4 (78,6 ^'8.7 
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'<- 




~«^i.o" 




K 




^2,0 




c% 




%0 




dl 




«4.0 




< 









fir 









^m 
^n 













_0 . 



(71) 



wIkto. 


~'^n{hb)yl^n{ha) 


ai.z= -j^^'nihc) 


^1.1=- 


-^niha) 




«1,2 = 


Uha) 


%.4=-^r.(M 


«1.3 = 


Mha) 


a,.= 'pUhO 


^2.0==" 


-l-^Ukib)^Uha) 


^5,7= p^A«(^3C) 


(^■2A = - 


-Y/n(ha) 


«'G.o= p ^niho) 


^2,2 = 


Y ^nihd) 


%6= Ylniha) 


^2.3 = 


Y/niha) 


^6.7= plA«(^l«) 


«;i.O^- 


-'Ukic)Mkib) ; 


«7.5=-^'An(*2a) 


«3.2 = 


Uhc) 


«7.6= ^ f^X^l«^) 


«^3.3 = 


Uhc) 


^7.7= ^^UM 


«^3.8= 


*r^"*"' 


%4=-Q.^.r«(M 


%o=- 


-^ rn{ha)^n{hb) 


a8.a= ^V^U/^ic) 


«4.2 = 


l-^ Cihc) 


(72) 


Under certain conditions, such as long wavelengths [Johler and Walters, 1960; Johler, 1961 ; I 
Jollier, Walters, and Harper, 1960], the effect of coupling in the region between r=a and r=c, 
figures 2, 3, and 4, can be neglected. This assumption reduces the matrix eq. (71) to, 




'-^niha) Liha) 


Mfcio) 


- 




"a„" 




--Ukib)4'n{kia) 




1 


-l^^'niha) l^jUfcia) 


j^ ^''nik^a) 






h„ 




^ ln(k,hWn{ha) 






Uhc) 


Mkia) [- 


-f«(^3C)] 




c„ 




- ^uikiC)Mkib} 






^ rU^i^) 


^m.c) 


-^ f;(^30)]_ 




d„ 




~r^ ^'n{kic)<P„ ik,b} 

_ "i J 




where a«=a^, bn=bn, Cn=c%, d, 
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(73) 



I 10-'^ 



* 






AMPLITUDE lErl 
NO IONOSPHERE 


^\:^ 


000 c/s 




€2- 

lo 


ID. Cr= U.UUD 

2 = 1 


: lOc/s 


Viooc^ 


k^ 




' ic/s : 


o^ 


L,.^^^^ 




io.c^- 


- 






X 


^ 100 c/s^ y 


: 




1 




1000 c/s 


lJ\ 



2000 



6000 
STATUTE MILES 



Figure 6. Amplitude of the vertical electric field, 
|Er|, of the ground wave excited by a vertical electric 
time harmonic dipole source current moment as a 
function of distance, d, along the surface of the 
terrestrial sphere, illustrating the standing wave at the 
antipode of the transmitter resulting from two waves 
traveling "around-the-world'' in opposite directions. 

/oZ=— --3.34(10-2) 
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PHASE LAG.-ARG E^ 
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NO IONOSPHERE 










«2 = I5, o- = 
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1000 


c/s^ , 


': 


-^ 
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- 
/I— — — ^ 










10 C/b 


iZJ^ : 
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Figure 7. 



Phase of the vertical electric field, arg Er, 
of the ground wave. 
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Figure 8. Amplitude of the horizontal magnetic 
field, |H,^|, of the ground wave excited by a vertical 
electric dipole time-harmonic dipole source current 
moment as a function of the distance, d, along the 
surface of the terrestrial sphere, illustrating the stand- 
ing wave near the antipode of the transmitter result- 
ing from two waves traveling '^around-the-world^' in. 
opposite directions. 

/iOC 



PHASE LAG, -ARG H^ 
;N0 IONOSPHERE 
i<2=l5, CT = 0005 
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Figure 9. Phase of the horizontal magnetic fields 
arg H^, of the ground wave. 



This matrix (eq. 73) can be generalized for straification of the sphere and the continuoulsy 
stratified concentric plasma, 



^1,1^1.2^1,3 

^2, 1^2, 2^2, 3 

^3, 2X3, 3^3, 4^3, 5 
^4, 2^4, 3^4, 4^4, 5 

^5, 4^5, 5^5, 



^6, 4^6, 5^6, 6^6, 7 

^7, 6^7, 7^7, 8^7, 9 
^8, 5^8, 7^8, 8^8, 9 



^2m-|-3, 2'W+2^2m + 3, 2m-f 3^2w + 3, 2^+4^2^ + 3, 2w + 5 
^2m+4, 2w-f 2^2w + 4, 2m-H3^2m-|-4, 2w+4^2m4-4, 2m + 5 



^2p4-3. 27? + 2^22? + 3, 2p+3^2p + 3,2p+4 
^2p-|-4, 2p + 2^2??-l-4. 2p + 3^2p + 4,2p+4_ 





rO.2, n 




~ "" 




An 








dhn 








bn 




^3,0 




Cn 




^4,0 




di,n 




^5,0 




ei,n 




^6,0 




d2,n 








dm-\,n 








^m — 1, n 








dm,n 








^m, n 








dp-i,n 








ep-l,n 








jlv,n J 




_ _ 



(74) 
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^here 



Xl,l=-^n(h,2f) 


X4.4=^ tniha) 


Xl,2=^n{h,J) 


X4,5=J^i^n{h(l) 


Xi,z=Mh,J) 


X5,0^-^n(hCl)Mf^lb) . 


X2.1 = -, 'An(^2.2/) 
/t2,2 


X5,4=fn(^lCl) 


K^2,\ 


X5,5=MhCl) 


X2,3 = l— i^Uhlf) 
/«^2,1 


X5,,= -^n(h,lCi) 


Xs,o=-^n{hb)Mha) 


X5.7=—i^n{h,lCl) 


X3,2=-fn(fe.l«^) 


x^,o=—j^ ^nihciyMhb) 



X3. 3 = — fn (h. iGi) Xs, 4=77^^ ihCi) 



X3. 4-- f n (^i«) Xe, 5=TrV'^ (^iCi) 






X3,5 = '^n(hCi) X6.6=— 77- ^n(h,lCl) 

X4,o=— 77 ^n{hb)rpn(ha) Xg.y^ — 77- i^nih.iCi) 

1^1 n^3. 1 

X4.2 = — ^ r^fel^^) X7.6=f;.fe,lC2) 

X4.3=— 7— ^n{h,l(j) X7,7 = ^^(A"3,lC2) 

'^2,1 

X7 ,8= — f ;i (^3. 2C2) 1_ ^/ / 7 X 

^3,m-l 
X7.9=—^/t (^3.2^2) 

1 X27/i+4,2m+3=T; 4^n\J^Z,m-lC'm) 

Xs Q=T fwW,lC2) 
"^3,1 

X2m+4.2w+4^= T fwW.T^^m) 

1 »r/7 \ ^3.m 

X8.7=T~^«(%1^2) 

"^3,1 1 

X2w+4,2m+5^= T Yn\'^Z,ni(^m) 

X8,8= T fnW,2C2) 
'<^3,2 

'^^'^"""'KT ^^^^3,2^2) X2p + 3.22, + 2=r;^(A-3,^,-lC^) 



3,2 



X2p+3,2p + 3 — Yn\'^Z,p-lCp) 
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^2p-|-3. 2p + i — r« (n^3, p^p) 



^2w + 3.2w + 






^2m+3,2w+3 — 'AnW. m-l^w) .. 1 ,/,7 n 

A2p + 4,2p + 3 T rnW,p-l('-p) 



^2w + 3, 



2W + 4 — T/iV^S, /w^w) 



— 1 , 



'^2m + 3,2m + 5 — YnK'^Z, m^m) f<^3,p 

The inner sphere, figure 4, has wave number k2,2, radius/, and Hertz vector, 



(75) 



n2.2=rr^S {2n+l)a2,n^Pnik2,2r)Pn (cose); 

ft2,2'0 n=0 

the outer layer of the earth has wave number k2,i, radius a, and Hertz vector, 
^'^'^==1^-^ ^ (2'^+l)K.V'.(^2.ir)+/i..rnfe.ir)]P, (cos d); 

i^2A'^ w=0 



f = 100 c/s 



REAL 

IMAGINARY 




Figure 10. Index of refraction, 7; = Re7; + iIm77, of 
the model plasma for various electron densities, N, 
and classical collisions, v, in the absence of a super- 
posed magnetic field, Hni = 0. 
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XGURE 11. Index of refraction, 77o = R'er7o + iIm7;o, 
of the model plasma for various electron densities, N, 
and classical collisions, v in the presence of a super- 
posed magnetic field, l^rn = 0.5 gauss, employing the 
quasi-longitudinal approximation. 
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Figure 12. Index of refraction, 77e = Re7;e+iImr;e, 
of the model plasma for various electron densities, N, 
and classical collisions, v in the presence of a super- 
posed magnetic field, 11^ = 0.5 gauss, employing the 
quasi-longitudinal approxiination. 
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the air (ki) is as before; the tnih layer of the ionosphere has wave number ^3,,,,, lower edge at 
r=Cm, <^iid Hertz vector, 

1 

Il3,m = r^ 1 S (2/^+1) [dn.m^n fe.m^O + ^n. m V'nC^S. m/') ]^^.. (cOS 6)] 

f^^.mf'O 11=0 

except for the outermost hxyer with wave number kz,p, which is infinite and thus has only 
outgoing waves, i.e.. 



^^'P==iJ-^ ^ ('^n+l)d,,, Uh,pr)Pn (cos e). 

1^3, p'O 71=0 



In each case, Il3= 1X3,0 + Il3,e. Similarly, eq (71) can be generalized to take account of coupling 
at each ionosphere boundary. 



(2n4i) 7^ 
f = 100 c/s 
h = 67.5 km 
N=56, 1/ =1.6x10^ 




^ -10 



rno 20 50 40 50 60 % 80 90 ioo Oo 120 jso WO iso 160 7 170 180 mm 





-100' — 

Figure 13. Locus of the coefficient eq {102), of the Legendre function Pn(cos Q) for spherical waves of order n mapped 
into the (2n-f-l) ja-plane, illustrating both amplitude and phase of the complex spherical wave coefficient. 
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2.5. The Watson Transformation 



The method proposed by Watson [1919] consisted of forming a complex integral for the 
Hertz vector, n, 



ni^T^^S (27l+l)P,(C0S^) MjL^^^^^+^^'n'Pnik.r)] 



D„ 



(76) 



thus, 






around a suitable contour, C, and taking the sum of the residues at the zeros, 

A-i=0, 



Ur- 



-27r^ y, -P%-i(-cose)[6;^_l^^_i(V)+<-j\A.,-i(V)] 



Arir6 ■^ cos j'^tt 



[l-M..,. 



+ 



f*ico 



IO-»rr 



10-8 



5 
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(79) 
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FiGUEE 14. Amplitude of the vertical electric field, ]Er|, at the surface of the terrestrial sphere (r = b = a) excited hy 
a vertical electric time harmonic dipole current source located on the surface of the terrestrial sphere at various 
frequencies, as a function of distance, d = a8, on the surface (r = a) of the sphere. 
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where the Hne integral along the imaginary axis is usually considered to be negligible, J{v) is 
defined by eq (77), and where 






^n — r% ' 



(81) 
(82) 



and Dn is the determinant of the matrix eq (73), for example. Whereas the series, Hiivs), con- 
verges rapidly in most practical computations, it would be necessary to evaluate the roots of 
the eq (78), and form the derivatives of the f^-waves and ^v-waves with respect to complex 
order, v. to maintain the equivalent rigor of the series of zonal harmonics. The exact equiva- 
lence of this procedure is of course assured by the theory of functions of the complex variable. 



1000 


^ 


1 1 

PHASE 

€ = 15, 


1 1 1 
LAG, -ARC 
a =0005 


1 1 


1 I 1 1 1 1 


_l_^^ 




z 


h 


675 


km 


^^,^-- 


■ ' 


- 




I 


N 


56. 


I/=l.6xl0\^^ 


'"^ 




I 




^ 






y^ 






- 


100 


1 


/ 


/ 




- 


^^__jooc/S--— 


~ 


10 




/ 


/ 






lOC/S 


—— 




z 


r 








IC/S 


I 




-i 
















-1 












- 






1 


1 


1 1 1 


1 1 


, 1 i 1 ^ 1 


1 1 1 1 I 



e" 


1 1 I 1 

PHASE LAG, 


1 1 1 1 

-ARG Er 


1 1 1 1 1 1 1 1 


\^J^-^^— 


- 


/^ 




\VX7S_ 


- 


1 


f ^ 


^ 


^ 


E 


L 
-I 


^ 




lOC/S 










IC/S 


- 


A 


ill! 


1 i 1 1 


1 1 1 i 1 ! 1 i 


1 1 1 



500 1000 1500 2000 2500 3000 3500 4000 4500 5000 
STATUTE MILES 

Figure 15. Phase of the vertical electric field, arg 
Er, and arg H</, at the surface of the terrestrial sphere 
as a function of distance. 



3. Computations 

In the particular case, figure 2, of a finitely conducting sphere (without a concentric 
plasma), the total Hertz vector in the air for the restriction, r=a, b>a, 



IS, 



ntot=no+nj, 
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(83) 
(84) 



where. 



7n = 






(85) 



where the Wronskian identity, [Watson, 1918], 



^n fn 



= — ^, 



(86) 



has been employed. From the definitions (62, 63) oi\l/n{^) and ^^^(2:) and the identity [Watson, 
1958], 



it is easily shown that 



j;(3)=J._l(2)-Jj.(2), 



^P'„{z) = ^u-l{z)--4'n{z). 



(87) 
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Figure 16. Amplitude of the horizontal magnetic field, |H^1, at the surface of the terrestrial sphere (r = b = 

a) at various frequencies as a function of distance. 
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Similarly, 



^n(^) = ^n-l(^)--L{^). 



(S8b) 



Thus, 7„ can easily be coinpiited by evaluating the functions rl/ni^) ^^^1 ^ni^)- Now, 

\l/_i{z) = cos z, 
and 

271+1 



Similarly, 



and 



fl(3) — 6--(l-i/2), 



2n+l 



^n(2) = L + l{z) + ^n-l(^)- 



(89) 



(90) 
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Figure 17. Am/plitude of the vertical electric field, |Er|, at the surface of the terrestrial sphere (r = b = a) excited\h'if 
a vertical electric time-harmonic dipole current source located on the surface of the terrestrial sphere, as a function 
of frequency, f, at various distances, illustrating an ^^ absorption band'' between 300 c/s to 3 kc/s. 
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Therefore it would appear that ^A«(^) and ^n(^) could easily be found for all n. However, for 
7i> 1^1, eq (89), loses precision with each recursion [Goldstein, and Thaler, 1959], and the series 
(84) does not converge until n:> \kia\. For 7i> \Icia\j the following method of Goldstein and 
Thaler [1959] is employed: 

For M> 1^1, define Fm+i(^)-0, i^M(^) = l, and i^n-i(^)=^^^Xi^n(s)-i^.+i(^) for n<M. 

3 



Then, with fixed L>l,n< (M-L), 



^„(3)=aF„(2)+6(M), 



(91) 



where a depends only on z and e can be made arbitrarily small by increasing M. a can be found 
from 






(92) 



for some A^ for which xpn is known, and then \l/n(^)y N<Cn<i{M-L), can be computed using eq (91). 
Since Im (^2^) is very large, yl/nikoo) and yf/' nQc^a) are very large. Therefore, it is more 
convenient to compute 
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Figure 18. Amplitude of the horizontal magnetic field, |H^| as a function of frequency for various distances, 
illustrating an "absorption band'^ between 300 c/s and 3 kc/s. 
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using, 



and, 



Rn = 



i?o~cot {k2a) 



n—k2aRn-i A"2« 



(94) 



(95) 



for n>0. This recursion formula, which is derived from the well-known recursion formulas 
for Bessel function [Watson, 1958] will lose precision for n^ \k2a\, but for the computations 
presented in this paper, n remained smaller than \k2a\. 

The Legendre polynomials P„(a:) can be computed from the recursion formula [Hobson, 
1931], 






(96) 



starting with P(o/) = l and Pi(x)=x. 

The convergence of the series (84) is notoriously slow, but the sum can be found from a 
more quickly convergent series. For ri^lkial [Watson, 1919], 



(2n+l)7. 2kia 



(97) 



ForZ<l, [Hobson, 1931] 
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PiGURE 20. Amplitude of the vertical electric field, |Er|, at 100 c/s as a function of distance from the source, illus- 
trating the primary field, the field disturbed by the terrestrial sphere {no ionosphere) , and the field disturbed by both 
the terrestrial sphere and the concentric plasma. 
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J2 /^PJcosi9)-(l-l/cos^+Z2)-^ 



so that 



n=Q 



k s ^^- (f)" ^"(-^ '^-kjb-' (f) ^«^ ^+9' • 



^a6 



(98) 
l99) 

(lOO) 



Therefore, adding (84) and (99), 

n-=A^ [^-T <-^ ^+9" -^ S [2^'^^ (f)%(2.+l)..] P„(cos 6). 
As a result of eq (97), 

as n gets much larger than kia, so that the new series (100) converges more quickly than the 
original series (84). The Er and H^ fields were computed at a given point, po, by computing 
the Hertz vector Iltot at seven points, j)_2, p_2, P-u P-o, Pu P2, P2, where Pt,=2h+^^A6 for 
fixed A^. The derivatives, eq (5), 

dd '^=^c 
and 
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PiGURE 21. Amplitude of the horizontal magnetic field, |H,^|, at 100 c/s as a function of distance from the source, 
illustrating the primary field, the field disturbed by the presence of the terrestrial sphere {no ionosphere), and the 
field disturbed by both the terrestrial sphere and the concentric plasma. 



/oZ=i^^3.34(10-2) 



763 



were then found numerically using the seven point Lagrangian formula [Kopal, 1955]. As a 
check, the first two derivatives of an orthogonal least-square polynomial fitted to points p_2y 
• ' ' y P2 [Hildebrand, 1956] were computed and compared with the results of the Lagrangian 
formula. The computed derivatives were then substituted into the appropriate eqs (5) to 
find the Er and H^ fields. 

For the model, figure 3, which includes the sharply bounded inosphere, computations 
were made for the case r=a=b, the Hertz vector for which can be written with the time 
function exp (ioot) suppressed: 



exp {—ikP) , 1 
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(101) 
(102) 



for the Gi^j defined in (72). 
Or, explicitly, 

yn=[(t>n[[hln{h(lWn{h(l)-hRnln{ha^^^ 

-(t>n { [kMnha) U (he) - klRn^Pn {ho) f ^ (^i^) ] [ f ^ (^ia)iA^ {kic) - U (he) ypn {ha) ] - ik2h U (he) i^n (ha) } ] 
X[0« { -k2[^n{ha)yPnhc) -^Pn(ha) Uhc)]+hRAUha)^Pn{hc) - ^Uhc)Mha)]} 
+<^n[hh[^'n(ha)yPn(he)-yl^n(ha)ln(hc)]-k\En^^ (103) 

where, 

<i>n=U(h,QC)+lri(k^^eC), 



f^3. n^3. e 



and, 

as before (93) . 



R 



1-3. e 
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The \l/ and f functions were eomputed as before, eqs (88-92). 

(V^rtaln nuinei-ical characteristics of the series (101) make possible direct computation 
of \\\v derivatives using 

"'"' ^ S (27^ + l)7.P; (COS d), (104) 



where PL (cos 6) denotes -7-7 — -— - P^ (cos d) . 
^ ^ d (cos d) ^ ^ 



Similarly, 



^=^2 S (2/^+1)7. [sin^ dP'^ (cos ^)-cos SP^, (cos 6)] 

X; (2n + l)y,Pn (cos 61)+-^-- -^ • (105) 



k'ia' t^o ' "" ' ' shi (9 d^ 

The derivatives of Pn(x) can be found from eq (96), 

PUi(.r)=P;_i(x) + (2/^+l)P.(^), (106) 

and [liobsoji, 1931] 

p':^,{x)=p':.,{x)+{27i+i)P',{x). (107) 

The series, eqs (104, 105), converge even more slowly than the original series, eq (101). 
However, for 7i>> \kia\, (2n+l)yn/^ikia, so, as before, the series can be transformed to 
improve convergence using [Hobson, 1931] 

- ri-P.(cos^) ^ 1.3....(2//^-l)A- 
,^n d (cos or (l-2Acos^+A2)^+^ ^ (/^<U Vi^»>> 

Substitution of tlie derivatives into eq (5) yields the Er and iJ,^ fields. 

The number of terms required for graphical accuracy in \Er\ varied from about 700 at 
10 c/s to about 3,000 at 10 kc/s. The CDC-1604 computer could compute about 830 7„, eq 
(102), per minute. Using these stored jn, the fields could be found for 40 different values of 6 
(distances) per minute per 1,000 terms required. 

At frequencies less than about 1,000 c/s, if low electron densities (<ClOO electrons/cm)^ 
are assumed, precision is lost during summation of series (105). This loss of precision increases 
with increasing d, and limits the usefulness of the zonal harmonics technique at ELF. 

4. Discussion 

A series of computations was carried out based on the model of a terrestrial sphere with 

47r 
a vertical electric dipole source current moment, eq (7), Io/= — ^ on the surface of the spliere.^ 

The field at the surface of tlie sphere, r=^a, figure 2, was described by the computation. Such 
a field is the exact equivalent of the ground wave described by van der Pol and Bienimer 
[1937, 1938, 1949] for example, utilizing the Watson transformation (section 2.5). Whereas, 
the works of van der Pol and Brennner [1937, 1938, 1949] employ a long series of approximations 



4/j=— ,^3.34(10"2) amp— ni. This is equivalent to the conventional notion of radiated power, Pr, as follows [Tohler, 1961]- Pr=l 6(10"3)a}2 
HoC • V / 

{loiyiZo, where Zo^VIOir (impedance of space). Such a radiation field cannot very readily be measured at ELF/VLF since the induction and 

.electrostatic fields are not negligible. 
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Figure 22. Amplitude of the vertical electric field, 
|Er|, as a function of distance for the terrestrial 
sphere with various model concentric plasmas at 
10 kc/s, (Hn, = 0, N = 0, 20, 56, 150) and 100 c/s 
Hm = 0, 0.5, N = 0, 0.25, 20, 56, 150) illustrating 
the effect of the earth's magnetic field and different 
electron densities. 
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Figure 23. Amplitude of the horizontal magnetic 
field, |H<^|, as a function of distance for the terrestrial 
sphere with various model concentric plasmas at 
10 kc/s, (Hn, = 0, N=0, 20, 56, 150) and 100 c/s 
(Rm = 0, 0.5, N = 0, 0.25, 20, 56, 150) illustrating 
the effect of the earth's magnetic field and different 
electron densities. 



to achieve numerical mastery of the very difficult formulas resulting from the Watson trans- 
formation, this solution does not emplo^^ any approximations, leaving only tlie assumed model. 
This solution is also equivalent to the ' 'acoustical shadow of a sphere '^ described by Lord 
Rayleigh [1904], and indeed some of the curves show a close resemblance. Rayleigh, however, 
did not treat the field at the surface of the sphere but at considerable distance. 

The regular decrement of the amplitude of the ground wave as a function of distance 
can be noted, figures 6 and 8, for both the \Er\ field and the l^^l field, at frequencies of 1, 10, 
100, 1000 c/s. However, near the antipode of the transmitter, cf= 10,000 to 12,500 statute 
miles, an interference pattern or standing wave can be noted as a result of the wave traveling 
around the earth in the opposite (-6) direction. The corresponding phase curves, arg Ery 
arg H^, are illustrated, figures 7 and 9. Whereas, the conductivity and dielectric constant 
are assumed to be finite, (e2==15, a-=0.005 mhos/meter), the second term in the denominator 
of eq (78), 



lC2\pn{hCi) 



^n{ha) 



« 



^nihd) 



is quite small compared with the first term (of the order of 10^ less) at 1 c/s and 10 c/s. The 

conductivity terms become significant at the frequencies greater than approximately 500 c/s. 

The introduction of a concentric plasma, figures 3 and 4, is quite intricate. However,. 

to gain some insight into the propagation of waves at great wavelengths, and demonstrate 



766 



I I I 



1~ 



I I ' I 



AMPLITUDE. |E,| 
€,M5, CT = 0.005 
H„=0.5 GAUSS. Io/ = 



MODEL [ h=67.5km, N = 56, j/= 1.6x10' 

MODEL 2 h = 75km, N = 160. i/ = 45xl0' 

MODEL 3 h = 80km, N = 500, i/ = I 6 x 10* 

MODEL 4 h = 90km, N = 10,000, ;/=L55xl0* 




1500 2000 2500 
STATUTE MILES 



3000 3500 



FiGUEE 24. Aniplitude of the vertical electric, |Er|, 
field as a function of distance, illustrating various 
reflection heights, h, of the concentric plasma, c = a 
+ h, and corresponding tj/pical ^^ daytime-noon^^ 
model plasmas {electron densities, N, and collision 
frequencies, v) . 
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Figure 25. Amplitude of the horizontal magnetic, 
IH^I, field as a function of distance, illustrating 
various reflection heights, h, of the concentric plasma, 
c = a + h, and corresponding typical "daytirne-noon^' 
model plasmas {electron densities, N, and collision 
frequencies, v), 
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the coinpiitational value of zonal liannonics, tlie particular case, figure 3, of a terrestrial sphere 
of finite conductivity surrounded by a concentric plasma was evaluated for the special situation 
in which transmitter and receiver were located on the surface of the sphere, hi=h2=0, as 
was illustrated for the ground wave case. The index of refraction of the concentric plasma 
was calculated, neglecting the earth^s magnetic fields, Hm=0, figure 10, and introducing the 
earth^s magnetic field, Hm=0.5 gauss with the aid of the Q-L approximation, figures 11 and 12, 
for various collision frequencies, p, and electron densities, A7'(electrons/cm^). Figure 13 illus- 
trates the behavior of the complex coefficient (2?i-f 1) jn, eq (102), of the Legendre polynomial 
P„(cos 6). Typically, the '^loops'' in the path occur for n^kia, the important terms are those 
for k^a^n<C2kia, and (2/1+1) 7^ approaches ik^a as n approaches infinity. 

The introduction of the concentric plasma causes a drastic change in the detail of the propa- 
gated field, figures 14 and 16. The space between the concentric plasma and the terrestrial 
sphere, ki, figure 3, is now filled with waves trapped by the two boundaries, r=a, r=c. Thus, 
figure 5, there are upgoing spherical waves, f-waves, and downgoing spherical waves, i/'-waves, 
whereas in the medium above the outer boundary, r^c, there are only upgoing, fo- waves and 
fg- waves (ordinary and extraordinar}^? ^^nd in the medium below the inner boundar}^ tliere are 
only ingoing i/'-waves. These spherical waves traveling in the radial direction, upon summa- 
tion, build up in the medium between r=a and r=c waves traveling the direction (±^) of in- 
creased angular distance. 

The regular nature of the field curves, \Er\ and \H^\, as a function of distance, d, figures 14 
and 16, indicates a single wave propagated in the ^-direction around the terrestrial sphere in 
the medium ki. Near the antipode of the transmitter, however, waves traveling in opposite 
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directions can be noted as a standing wave, with a half wavelength between minima, as was 
noted in the case of the ground wave. At the lower frequencies, 100 c/s, and higher conducting 
concentric model plasmas, figures 24 and 25, this phenomenon can be noted at short distances 
«4,000 statute miles). Similar standing wave phenomena can be noted at the higher fre- 
quencies, 10 kc/s, figures 24 and 25, and 30 kc/s, figures 14 and 16, for intermediate conductivity 
concentric model plasma, indicating interaction between two or more waves. The distance 
between minima are in this case much greater than a wavelength, and seem to damp out with 
distance. 

The variation in the amplitude, |E;.|, |i?</,|, figures 17 and 18 and the phase, figure 19, as a 
function of frequency is illustrated. These curves illustrate the spectrum mutilation [Johler, 
1962] of the pulses propagated through the media defined by this paper. It is interesting to 
note an ''absorption band," between frequencies 300 c/s and 3,000 c/s, for the assumed model, 
figures 17 and 18. This absorption band is also reflected in the phase correction curves, figure 
19. The primary field phase is also illustrated. It should be noted that phase is presented as a 
lag, — arg Er, — arg i?^, or — arg Er/Er, pri, — arg H^/H^,, pri, where Er, pri and H^^ j,ri refer to 
the primary fields. 

A comparison (at 100 c/s) of the primary field, the ground wave (no ionosphere), and the 
field with the concentric plasma (ionosphere) is illustrated for the E'^.-field, figure 20, and the 
7^^-field, figure 21. For the assumed model, the attenuation rate of the field increases succes- 
sively as the terrestrial sphere and the concentric plasma are introduced. 

It is evident, figures 22 and 23, ^"^^^^^le electron density, N, is the most important con- 
centric plasma (ionosphere) factor governing the field. The intensity of the superposed mag- 
netic field (earth's magnetic field) is only secondary. It should be noted, figures 22 and 23, that 
the field reduces to the ground wave for N=0. The vertical electric field lE^I and the horizontal 
magnetic field \Hcf>\ are illustrated. Model ionospheres employed in other analyses [Johler and 
Harper, 1962] by selecting the electron concentration values of the earth's ionosphere at day- 
time noon, together with appropriate collision frequencies, are illustrated, figures 24 and 25. 
It is interesting to note the increase in the field at 100 c/s as the higher conductivity models are 
successively employed (model 1, 2, 3, 4). However, an intermediate model 2 exhibits not only 
the greater field at 10 kc/s but also the standing wave with distances between minima much 
greater than a wavelength. 

Since, analysis is designed to gain some insight into the very complicated behavior of 
terrestrial radio waves near the earth, and also to demonstrate the technique of zonal har- 
monics as applied to the radio problem employing modern techniques of computation, con- 
siderable caution must be exercised in the interpretation of results from a practical viewpoint. 
The roughness, daytime-nighttime changes, the temporal variations and the multiplicity of 
frequencies in a radio signal, and coupling at the ionosphere all could tend to mask the theo- 
retical phenomena described. Furthermore, it would be desirable to investigate the effect of 
the heavy ions numerically. But the detailed investigation of these matters is reserved for 
future studies. 

5. Conclusions 

The results of this paper demonstrate that the field of the low frequency or long wavelength 
radio wave, less than approximately 50 kc/s, can indeed by evaluated by a summation of a 
series of zonal harmonics. Whereas the number of terms in the series could become quite 
large, the speed with which these terms can be summed on a large scale computer offsets the 
complications introduced by the Watson transformation. 

The detailed structure of the field is characterized in the absence of a concentric plasma 
(no ionosphere) by the quite regular behavior of the ground wave as a function of distance. 
Indeed, the steady decrement of the ground wave field is modified only near the antipode, 
whereupon an interference pattern or standing wave as a function of distance is noted as a 
result of interference between waves traveling ^'around the world' ^ in opposite directions from 
the transmitter. 
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The introduction of the concentric electron-ion plasma shell traps the waves leaking 
into space, which upon reflection from the plasma build up traveling waves in the direction of 
increased distance from the transmitter. Thus, the series of zonal harmonics comprises 
individual waves which are traveling in the radial direction with respect to the center of the 
sphere and standing in tlie (hrection of increased anguhir distance around the sphere. These 
waves when sununed build up the wave progressing in the direction of increased angular dis- 
tance. But under certain special conditions such splierical waves also build up a standing 
wave in the same direction which can be interpreted as a multiplicity of traveling waves such 
as two waves traveling in opposite directions at the antipode of the transmitter. 

The results of the computations indicate that full rigor can be achieved with comparative 
ease at frequencies less than 50 kc/s, leaving only the assumed model for the transmitter and 
the propagation medium and avoiding the complications of the Watson transformation. 

The results of this paper suggest refinements of the computations, such as the introduction 
of terrestrial sphere stratification, continuous ionosphere stratification, polarization coupling at 
the ionosphere, and the effect of heavy ions. 
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7. Appendix 

The eoefRcients of the quartic (23) are: 

/ 2 1x2/1 '^' \ I / 2 i\ /l 1 ''^— 2 aLhr^ \ . s—1 
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where , 






h= — ^; 



hL=--—h sin /, a=sin <^i, 

hr^h cos /, co%=Ne^/Km, 
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ax=sin (t)i cos 0„, coH^fioeHJrn, 

a7'=sin 0/ sin 0„, K==l/c^'/xo=eo. 

These coefficients can Ix* modified lor (4ectroii-ion collisions proportional to (Mi(»r<i:y ['Tohlcr 
and Harper, 1962]. The angle of incidence, eq (47), (/), = sin~^ ^ h ^"^^^ ^^^' found i'oi' (^ich (^1(»- 
mentary spliorical wave, thus, the* f,rwave can })e wj'itten for two wave types, 1, 2, 



U.J .,2 

and the i/^,rwave, 
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(t>i=sin '— f 



Ke 



a complox an^lo of incidcnco, or, 



for either \j/- or f-waves and (^ is the second solution of eq (60). 

The Q-L appi'oxiniation can be replaced by the full magneto-ionic theory by evaluating the 
index of refraction, 77, Irom eqs (23 and 24) and employing th(» exact expression for P and Q in 
eqs (49, 50, 51, 52, 53) [Johler and Harper, 1962]; 
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The coefficients of the matrices in eqs 31 and 32 are: 
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(i2z=—(iz2= __ y =Q? if Hm is in the 2/^-plane 



W+ /^ ^4^12\ 

m+ / . ^ C4(7n\ 
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(246^ ae4= 
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cf56=— 065= — \r^=Oi if ^m IS in the yz-plane 
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%9 — ^97 = 
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^89=— «98=3-i^=0, if H,n is in the yz plane 
where the C* are defined in eq (29) and Si=Cj+Ci+i + C'i+2. 

^12= — OjifLji 
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b2i=~aT(iL 

^23=— fl^'Lf 
632 = 05-^0' T' 

633=ay2 + araL— 1. 

It is clear that a fundamental difficulty exists in treating anisotropic media with the wave 
eq (55). Thus, the separation of variables, eq (56), depends upon the assumption of isotropic 
media of propagation, i.e., /:3 must be independent of position and direction of propagation. 
Indeed, the complications introduced by such considerations have only been satisfactorily 
resolved in plane coordinates b}^ treating Maxwell's equations directly. The procedure which 
has been outlined in this paper introduces such solutions of the nuigneto-ionic theory ad hoc 
to the rigorous series of zonal harmonics. Tliis is a specification for the plasma medium k^ 
which is an approximation. The validity of this approximation at long wavelengths depends 
on the behaviour of the complex index of refraction with respect to direction of propagation, 
such as the ratio, Ee,n/Er,n, for elementary spherical wave, which can be readily computed 
with the aid of methods described in this paper. Also, the index of refraction, for sin ^i= 
sin ^i^n dependence on and n can be evaluated. If n <^kiC, Ee,n/Er,n for each wave type 1, 2 
is almost real. Also, if n^kiC an elliptical polarization can be noted, since [Ee^n/Er^nlia are 
each complex. Furthermore, if ri^kjC, each [Ee^n/Erjilio becomes imagiiiaiy indicating circular 
state of polarization. 

The introduction of these complex numbers, sin cj^i, into the quartic eq (23) results in an 
approximately constant index of refraction, rj for electron densities, A^=50 — 1000 electrons/ 
cm^ over all distances, d^ad, for values of n^O to 10 kiC, ytiC-^ 10 — 1000. This is especially 
true in the more important region n^kiC. This is consistant with the demonstrated validity 
of the Q—L approximation at long wavelengths [Johler and Walters, I960]. The degree to 
which the method of introducing the magneto-ionic theory to the rigorous series of zonal har- 
monics can be extrapolated to shorter wavelengths or higher frequencies is worthy of further 
investigation. Otherwise, it is necessary to return to Maxwell's eq (1) and reformulate thc^ 
problem. It is anticipated, however, that the essential features of the magneto-ionic theory, 
such as the effect of the direction of the earth's magnetic field vector, the effect of polarization 
coupling external to the ionosphere, and the effect of the heavy ions can be adequately treated 
at long wavelengths with this theory. 

(Paper 66D-234) 
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